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We propose a conjecture for the exact expression of the dynamical zeta function for a family 
of birational transformations of two variables, depending on two parameters. This conjectured 
function is a simple rational expression with integer coefficients. This yields an algebraic value 
for the topological entropy. Furthermore the generating function for the Arnold complexity is also 
^ • conjectured to be a rational expression with integer coefficients with the same singularities as for 

the dynamical zeta function. This leads, at least in this example, to an equality between the Arnold 
complexity and the exponential of the topological entropy. We also give a semi-numerical method 
■ to effectively compute the Arnold complexity. 

O PACS numbers: 05.45. +b, 47.52. +j 
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I. INTRODUCTION 

To study the complexity of continuous, or discrete, dynamical systems, a large number of concepts have been 
introduced |§]|§]. A non exhaustive list includes the Kolmogorov-Sinai metric entropy the Adler-Konheim- 

McAndrew topological entropy the Arnold complexity Q, the Lyapounov characteristic exponents, the various 
fractal dimensions, [0J|] the Feigenbaum's numbers of period-doubling cascades |5|,[l0|, ■ ■ ■ ■ Many authors have tried 



to study and discuss the relations between these various notions in an abstract framework [11 12 1. Inequalities have 
been shown, for instance the metric entropy is bounded by the topological entropy, let us also mention the Kaplan- 
Yorke relation |l^,Q. Furthermore, many specific dynamical systems have been introduced enabling to see these 
notions at work. Some of the most popular are the Lorentz system |1| , the baker map [^6) , the logistic map jL7) , 
the Henon map |l8[ . Each of these systems has been useful to understand and exemplify the previous complexity 
measures. 

Here, we introduce another two-parameter family of mapping of two variables, originating from lattice statistical 
mechanics for which much can be said. In particular, we will conjecture an exact algebraic value for the exponential of 
the topological entropy and Arnold complexity. Furthermore, these two measures of complexity are found to be equal 
for all the values of the two parameters, generic or not (the notion of genericity is explained below). A fundamental 
distinction must be made between the previously mentioned complexity measures according to their invariance under 
certain classes of transformations. One should distinguish, at least, two different sets of complexity measures, the ones 
which are invariant under the larger classes of variables transformations, like the topological entropy or the Arnold 



*e-mail : dauriac@crtbt.polycnrs-gre.fr 
t e-mail : hassani@ist.cerist.dz 
*e-mail : maillard@lpthe.jussieu.fr 



1 



complexity 0], and the other measures of complexity which also have invariance properties, but under a "less larger" 
set of transformations, and are therefore more sensitive to the details of the mapping (for instance they will depend 
on the metric). 

We now introduce the following two parameters family of birational transformations k a y. 

u n +i = l — u n + u n /v n 

v n+ i = e + v n - v n /u n + a ■ (1 - u n + u n /v n ) (1) 
which can also be written projectively : 

^n+l iPn^n ^nVn ~t~ ^n^n) ' 

v n +\ = e • u n ■ v n ■ t n + (u n - t n ) ■ v n + a- (v n t n - u n v n + u n t n ) ■ u n 

tn+l = U n -V n -t n (2) 

As far as complexity calculations are concerned, the a = case is singled out. In that case, it is convenient to use a 
change of variables (see Appendix |A]) to get the very simple form k e : 

Vn+l = z n + l-e 

Zn+l = Vn ■ ^ ,\ (3) 

or on its homogeneous counterpart : 

yn+l ~ ~~ £ t-n) ' {%n ^n) ) 

t n+ i = t n ■ (z n + t n ) (4) 

These transformations derive from a transformation acting on a q x q matrices M [p(if : 

K q =toI (5) 

where t permutes the entries Mx$ with M^^, and I is the homogeneous inverse: I(M) — det(M) • M -1 . Transforma- 
tions of this type, generated by the composition of permutations of the entries and matrix inverse, naturally emerge 
in the analysis of lattice statistical mechanics symmetries . 

II. THE COMPLEXITY GROWTH 

The correspondence |20| between transformations K q and k a . e , more specifically between K q and k a ,e, is given in 
Appendix [A|. It will be shown below that, beyond this correspondence, K q and k a ^ share properties concerning the 
complexity. Transformation K q is homogeneous and of degree (q — 1) in the q 2 homogeneous entries. When performing 
the n th iterate one expects a growth of the degree of each entries as (q — l) n . It turns out that, at each step of the 
iteration, some factorization of all the entries occurs. The common factor can be factorized out in each entry leading 
to a reduced matrix M n , which is taken as the representent of the n th iterate in the projective space. Due to these 
factorizations the growth of the calculation is not (q — 1)™ but rather A™ where generically A is the largest root of 
1 + A 2 — A 3 = (i.e. 1.46557123 < q—1) |pO| , pl| . We call A the complexity growth or simply the complexity. This result 
is a consequence of a stable factorization schem^ given in Appendix |b], from which two generating functions a(x) 
and (3(x) can be constructed. Generating function a{x) keeps track respectively of the degrees of the determinants of 
the successive reduced matrices and (3{x) of the degrees of the successive common factors. The function a(x) should 
not be confused with the parameter a. The actual value of A is the inverse of the pole of j3{x) (or ot{xj) of smallest 
modulus. The algebraicity of the complexity is, in fact, a straight consequence of the rationality of functions a(x) 
and f3(x) with integer coefficients p0|. The same calculations have also been performed on transformations Eq. (PP) 



1 Complexity growth can also be understood from a singularity point of view ^7|, or through recurrence relations associated 
with the geometry of the singularities of the evolution [W8|. This is not the approach developped here. 
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and Eq. (g). In that case, factorizations also occur at each step, and generating functions can be calculated. These 
generating functions are, of course, different from the generating functions for if 2 (see f2fj|| ) but they have the same 
poles, and consequently the same complexity growth. One sees that, remarkably, the complexity A does not depend 
on the birational representation considered: for any value of q, k a ^ or the homogeneous transformation Eq. (|^). 
It will be usefull to define some degree generating functions G(x) : 

G(x) = ^ d„ • x " ( 6 ) 

n 

where d n is the degree of some quantities we look at, at each iteration step (numerators or denominators of the two 
components of k n , degree of the entries of the "reduced" matrices M n 's, degree of the extracted polynomials f n 's 
in Appendix [b] ...). The complexity growth A is the inverse of the pole of smallest modulus of any of these degree 
generating functions G(x) : 

lo gA= Um (7) 



A. Complexity growth for a = 

In the a — case, which corresponds to a codimension one variety of the parameter space (see Appendix ^ 
and ^|), additional factorizations occur reducing further the growth of the complexity. The generating functions are 
modified and the new complexity is given, for K q , by the equation: 

1 - A 2 - A 4 = (8) 

i.e. A ~ 1.27202 ■ ■ • . For k e , which corresponds to if 2 , the equation reads: 

1 - A - A 2 = (9) 

leading to the complexity A ~ 1.61803- •• ~ (1.27202 •• -) 2 . Not surprisingly, the complexity of the mappings k a ^ e 
for a = (see (|l|)) and mapping k e (see (||)), are the same: complexity A corresponds to the asymptotic behavior 
of the degree of the successive quantities encountered in the iteration (sec (vn). Clearly, this behavior remains 
unchanged under simple changes of variables. Note that this complexity growthn analysis can be performed directly 
on transformation k e , or on its homogeneous counterpart Eq.( ^). The number of generating functions in the two 
cases is not the same, but all these functions lead to the same complexity. In fact, complexity A is nothing but the 
Arnold complexity || , known to be invariant under transformations corresponding to a change of variables (like the 
change of variables from Eq. ([!]) (for a = 0) to Eq. (||) or to Eq. (0)). Let us also recall that the Arnold complexity 
counts the number of intersection between a fixed linefj and its n iterate, which clearly goes as A™. Conversely, all 
these calculations can be seen as a handy way of calculating the Arnold complexity. 

All these considerations allow us to design a semi-numerical method to get the value of the complexity growth A 
for any value of the parameter e. The idea is to iterate, with (^J) (or ([!])), a generic rational initial point (yo, Zq) and 
to follow the magnitude of the successive numerators and denominators. During the first few steps, some accidental 
simplifications may occur, but, after this transient regime, the integer denominators (for instance) grow like A™ where 
n is the number of iterations. Typically, a best fit of the logarithm of the numerator as a linear function of n, between 
n = 10 and n — 20, gives the value of A within an accuracy of 0.1%. An integrable mapping corresponds to a 
polynomial growth of the calculations : the value of the complexity A has to be numerically very close to 1 . Fig. ^ 
shows the values of the complexity as a function of the parameter e. The calculations have been performed using an 
infinite- precision C-library po]] . 

For most of the values of e we have found A ~ 1.618, in excellent agreement with the value predicted in Eq. (||). 
In (27|, it has been shown that the simple rational values e = —1, 0, 1/3, 1/2, 1 yield integrable mappings. For these 



2 Growth of the calculations related with factorizations were also introduced by Veselov for some particular Cremona trans- 
formations (^Q. 

3 Or the intersection of the n-th iterate of any fixed algebraic curve together with any other possibly different but fixed algebraic 
curve. 
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special values one gets A « 1 corresponding to a polynomial growth |^7j . In addition, Fig. [j] singles out two sets of values 
{1/4, 1/5, 1/6, • ■ • , 1/13} and {3/5, 2/3 , 5/7}, suggesting two infinite sequences e = 1/n and e = (m — l)/(m + 3)^| 
for n and m integers such that n > 4 and m > 7 and m odd. We call "specific" the values of e of one of the two 
forms above (together with the integrable values), and "generic" the others. To confirm this conjecture, we go back 
to (the matrix) transformation K q , for q = 3, to get a generating function of the degree s of some factors (the /„'s in 
Appendix [§|) extracted at each step of iteration, namely, with the notations of |4^j2l|]l9|] and of Appendix |§[ function 
(3{x). From now on, we will give below, instead of (3{x), the expression of the following complexity generating function 
defined, for q x q matrices, as : 

G°(q,x) = ^ (10) 
q ■ x 

In the following the calculations are often displayed for 3x3 matrices and Gf(q, x) will simply be denoted G" (x) . 
Let us recall that the value of the complexity A is the inverse of the root of smallest modulus of the denominator 
of this rational function. Examples of these calculations in order to get the corresponding factorization scheme and 
deduce the generating function (3{x) or G"(x) , are given in Appendix |b| In Appendix we show how to choose an 
initial matrix to iterate, the matrix satisfying a = and e = p/q for any integers p and q. First, we have obtained 
(see Appendix [b]) the generating function G e (x) in the generic case for a = : 

G e (x) = ] + \ +X \ (11) 
1 — ar — x 

We also got the generating function G c (x) for the different "specific" cases : 

„ . . 1 + x + x 3 - x 2m+1 - x 2m+3 . , „ 

G i/m(x) = 1 _ x2 _ xi + x2m+4 , with m>4 (12) 



G (m _ 1)/{m+3) (x) = 1 _ J .2_ X 4 + X 2m+A ' Wlth 171 ^ 7 171 ° AA ( 13 ) 

and : 



int[X) ~ l-x 2 -x 6 + x s -x w + X 12 + x 16 - x ls [ ' 

1 + X ■ (1 + X 2 ) + x A ■ (1 + x A + X 8 ) 



1 - x 2 ■ (1 - x 12 ) - x 6 • (1 - X 2 + x 4 - x & + x s - X 10 + X 12 ) 

for the two integrable values e — 1/2 and e = 1/3. For e = l/m (m > 4) and e — (m — l)/(m + 3) (m > 7 and m 
odd) , the corresponding complexities are the inverse of the roots of smallest modulus of polynomial : 

1-x 2 -x A -x 2m+i = (15) 

in agreement with the values of Fig. [l]. In this figure the e-axis has been discretized as M/720 (M integer) and the 
extra values 1/7, 1/11, 1/13 and 5/7 have been added. This semi-numerical method acts as an 'integrability detector' 
and, further, provides a simple and efficient way to determine the complexity of an algebraic mapping. Applied to 
mappings Eq. (|l|), Eq. (|^), or Eq. (||), it shows that the complexity is, generically, independent of the value of the 
parameter e, except for the four integrable points, and for two denombrable sets of points^. 

It is worth noticing that these results are not specific to 3 x 3 matrices, for example relation Eq. ([ll]) is actually 
valid simply replacing G"(x) by Gf(q, x) . 



4 Note that m — > (m + 3) /(m — 1) is an involution. 

5 These two sets of points also appear naturally in the framework of a "singularity confinement analysis" 
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FIG. 1. Complexity A as a function of e taken of the form M /720 plus the special values 1/7, 1/11, 1/13 and 5/7 for a = 0. 
The arrow indicates the expected value. 
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FIG. 2. Phase portrait of k e for a = and e = 13/25. 550 orbits of length 1000 have been generated. 50 orbits start from 
points randomly choosen near a fixed point of order 5 of k f = £13/25, and 500 others orbits start from randomly choosen points 
outside the elliptic region. Only the points inside the frame are shown. 



B. Complexity growth for a/0 



These complexity growth calculations can straightforwardly be generalized to a ^ 0. As explained in Appendix |B|, 
the "generic" generating function is : 
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1 - x 



(16) 



The pole of smallest modulus of Eq. (|16|) gives 1.46557 • • • for the value of the complexity for the matrix transformation 
K. The complexity for the transformation k a t is the square of this value: A = 2.14790- •• . Fig. || shows, for 
a = 1/100, complexity A as a function of the parameter e, obtained with the semi- numerical method previously 
explained. Even with such a "small value" of a the expected drastic change of value of the complexity (namely 
1.61803 — > 2.14790) is non-ambiguously seen. Moreover, Fig. || clearly shows that, besides the value e = known 
to be integrable whatever a |^] , at least the following values e = 1/2, e = 1/3 and e = 3/5 are associated with a 
significantly smaller complexity, at least for the discretization in e we have investigated. From these numerical results 
and by analogy with a = 0, one could figure out that all the e = 1/m are also non-generic values of e. In fact a 
factorization scheme analysis like the one depicted in Appendix |b|) shows that e = 1/4 ore = 1/7 actually correspond 
to the generic generating function Eq. (|l6|). We got similar results for other values of a ^ 0. However, when varying 
a and keeping e fixed, new values of the complexity A occur, A being some "stair-case" function of a. We will not 
exhaustively describe the rather involved "stratified" space in the (a , e) plane, corresponding to the various "non 
generic" complexities. Let us just keep in mind that, besides e = and e = — 1, at least e = 1/2, e = 1/3 and 
e = 3/5 are singled out for a =/= in our semi- numerical analyzis. The generic expression (for 3x3 matrices) for the 
generating function G{x) Eq. Jl6| ) is replaced, for the "non-generic" e = 1/2 (with a ^ 0), by : 
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For the other "non-generic" value of e, e = 1/3, the complexity generating function reads : 

r * M = l + x + x* - x 12 

^ ' (l-x 2 )-{l-x 2 -x 4 -2x G -x 8 -x 1 °) 
For the "non-generic" value e = 3/5, the complexity generating function reads : 

l + x + x 3 - x 20 



G V 5 ^ - (l-x 2 )-(l-x 2 ~ 



x 4 - 2x e -x s - 2x w - x 12 - 2x 1A - x 16 - a; 18 ) 
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FIG. 3. Complexity A as a function of e taken of the form M/720 for a — 1/100. The arrow indicates the expected value. 







III. DYNAMICAL ZETA FUNCTION AND TOPOLOGICAL ENTROPY 



It is well known that the fixed points of the successive powers of a mapping are extremely important in order to 
understand the complexity of the phase space. A lot of work has been devoted to study these fixed points (elliptic 
or saddle fixed points, attractors, basin of attraction, etc), and to analyse related concepts (stable and unstable 
manifolds, homoclinic points, etc). We will here follow another point of view and study the generating function of the 
number of fixed points. By analogy with the Riemann £ function, Artin and Mazur |2q| introduced a powerful object 
the so-called dynamical zeta function: 



( 00 fm\ 

C(t) = exp PT #fix(fc m ) ■ — (20) 

\m=l m J 

where #fix(fc m ) denotes the number of fixed points of k m . The generating functions 

H(t) = J2#Mk m )-t m (21) 

can be deduced from the £ function: 

H(t) = il(bgC(t)). (22) 
The topological entropy h is related to the singularity of the dynamical C function: 

lo gh = lim ^Wfc")) (23) 

rn— *oc 777 

If the dynamical zeta function can be interpreted as the ratio of two characteristic polynomials of two linear operators^] 
A and B, namely ((t) = det(l — t ■ B) / det(l — t-A), then the number of fixed points #fix(fc m ) can be expressed from 
Tt(j4") — Tr(B n ). In this linear operators framework, the rationality of the £ function, and therefore the algebraicity 
of the exponential of the topological entropy, amounts to having a finite dimensional representation of the linear 
operators A and B. In the case of a rational £ function, the exponential of the topological entropy is the inverse of 
the pole of smallest modulus. Since the number of invariant points remains unchanged under topological conjugaison 
(see Smale |}3| for this notion), the £ function is also a topologically invariant function, invariant under a large set 
of transformations, and does not depend on a specific choice of variables. Such invariances were also noticed for the 
complexity growth A. It is then tempting to make a connection between the rationality of the complexity generating 
function previously given, and a possible rationality of the dynamical £ function. We will also compare the Arnold 
complexity A and the exponential of the topological entropy h. 

A. Dynamical zeta function for a = 0, e generic 

We try here to get the expansion of the dynamical zeta function of the mapping k e (see Eq. ([}])), for generic 
values of e which are neither of the form 1/m, nor of the form (m — l)/(m + 3). We concentrate on the value 



e = 13/25 = 0.52. This value is close to the value 1/2 where the mapping is integrable |27|. On can gain an idea 
of the number, and localization, of the (real) fixed points looking at the phase portrait of Fig. 0. The elliptic fixed 
point (yo,zo) = (.24, —.24) is well seen, as well as the five elliptic points and the five saddle points of k^. Many 
points of higher degree are also seen. Transformation k e has a single fixed point for any e. This fixed point is elliptic 
for e > and localized at (yo, zq) = ((1 — e)/2, (e — l)/2). Transformation k\ has only the fixed point inherited from 
k e . The new fixed points of fcjf are (2 — e , (e — l)/2), (—1, 1) and ((1 — e)/2, e — 2). Transformation k\ has four new 
fixed points. At this point the calculation are a bit too large to be carried out with a literal e, and we particularize 
e = 13/25. For k\ we have five new elliptic points and five new saddles points. The coordinates z and y of these 
points are roots of the two polynomials: 



6 For more details on these Perron-Frobenius, or Ruelle-Araki transfer operators, and other shifts on Markov's partition in a 
symbolic dynamics framework, see for instance |29|- |3l| , fIof . 
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P(z) = (4375 z 2 + 1550 z - 89) (175 z 2 + 106 z + 7) (25 z 2 + 12 z + l) 2 (25 z + 6) 3 (24) 
Q(y) = P(-y) (25) 



The five pairings of the seven roots of Eq. (24) and Eq. (E5J), giving the five elliptic points, are (0.530283, -0.107335), 
(-0.050283, -0.24), (0.372665, -0.372665), (0.107335,-0.530283), (0.24, 0.050283) and the five pairings giving the five 
hyperbolic-saddle points are (0.372665, -0.075431), (0.107335, -0.107335), (0.404568, -0.24), (0.075431, -0.372665,), 
(0.24, -0.404568). This is clearly seen on Fig. || where the occurrence of five "petals" corresponding to five elliptic 
points are obvious, the five hyperbolic points being located between the petals. 

For transformation k® , beyond the fixed points of k and k 3 , one gets two complex saddle fixed points, i.e. transfor- 
mation k has two 6-cycles. For transformation kj , one obtains one elliptic real fixed point, one saddle real fixed point 
and and two complex saddle fixed points. For transformation kf , one obtains one saddle real fixed point and four 
complex saddle fixed points. For transformation k^ , one obtains one elliptic real fixed point, three saddle real fixed 
points and and four complex saddle fixed points. For transformation k\® , one obtains one elliptic real fixed point, 
one saddle real fixed point and and three complex elliptic fixed points and six saddle complex fixed points. The two 
elliptic fixed points of k\° (0.24, -0.874) and (0.874, -0.24) are seen as "ellipse" on Fig. (||). For transformation fc* 1 , 
one obtains one elliptic real fixed point, five saddle real fixed point and and twelve complex saddle fixed points. On 



Fig. (|2h a fixed poin t of k 12 lying on y + z = is seen near y — —13/25. The polynomials, similar to Eq. (24) and 
Eq. ( p5| ) (or to Eq. (Dl) given in appendix [E5|) , as well as the specific pairing of roots, for the successive iterates k N , 



are available in 

It is worth noticing, that among the 53 cycles of k e of length smaller, or equal, to 11, as much as 44 have a 
representent on the line y + z = 0, six have one on the line y + z = 0. Two of the three remaining cycles are of length 
11, while the last is of length eight. The particular role played by the y + z = line can be simply understood. Let 
us calculate the inverse of the birational transformation (|3j) . It has a very simple form : 

z n +i = y n - (1 - e) , y n+ i = z n ■ Vn +E (26) 

Vn - 1 

which is nothing but transformation (|^) where y n and — z n have been permuted. The ?/«<->— z n symmetry just 
corresponds to the time-reversal symmetry k e <-> k~ l transformation. The y + z = line is the time-reversal 
invariant line. Also note that only one of the 31 complex cycles is of the form Zq, Z\, ■ ■ ■ Z p , Zq 7 Z\, ■ ■ ■ Z p where 
%i = {l/ii z i) an( i %i is t ne complex conjugate. The 30 remaining complex cycles are actually 15 cycles and their 
complex conjugates. 

Eventually, we observe an area preserving ]35[ | property in the neighborhood of all the fixed points of fc" : the 
product of the modulus of the two eigenvalues of the Jacobian (i.e. the determinant) of fc™, at all fixed points for 
n < 11, is equal to 1 . This local property is rather non trivial : the determinant of the product of the jacobian over 
an incomplete cycle is very complicated and only when one multiply by the last jacobian does the product of the 
determinants shrink to 1 . 

The total number of fixed points of fc^ for N running from 1 to 11, yields the following expansion, up to order 
eleven, for the generating function H{t) of the number of fixed points: 

H t (t) = t + t 2 +4t 3 + 5t 4 + 11 + 16 1 6 + 29 i 7 + 45 1 8 + 76 i 9 + 121 1 10 + 199 i 11 +••• (27) 

This expansion coincides with the one of the rational function : 

H ' {t) (i-t^(i-t-t^) (28) 

which corresponds to a very simple rational expression for the dynamical zeta function: 

CM - (29) 



Expansion (27) remains unchanged for all the other generic values of e we have also studied. 

We conjecture that : The simple rational expression Eq. (|29| ) is the actual expression of the dynamical zeta function 
for any generic value of e. 

Comparing the expression Eq. (|J) with Eq. ((H), one sees that the singularities of the dynamical zeta function 
happen to coincide with the singularities of the generating functions of the Arnold complexity. In particular, the 
complexity growth A and the exponential of the topological entropy h are equal. 

When mentioning zeta functions, it is tempting to seek for simple functional relations relating £(t) and £(l/i). Let 
us introduce the following "avatar" of the dynamical zeta function : 
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(30) 



The transformation z — > z/(z — 1) is an involution. One immediately verifies that Ce(*) corresponding to (|29| ) 
verifies two extremely simple and remarkable functional relations : 

= -Ce(lA), and: £(*) = £("!/*), (31) 

or on the zeta function £(t) : 

Ce(1/<) = 2 -off- r and : u ~ l/t) = Ut) (32) 

The generating function (||^) verifies : 

H e (-l/t) = -H&) (33) 

An alternative way of writing the dynamical zeta functions relies on the decomposition of the fixed points into cycles 
which corresponds to the Weyl conjectures p6| ]. Let us introduce N r the number of irreducible cycles of k r e : for 
instance for Nx% we count the number of fixed points of k^epsilon 12 , that are not fixed points of k e , k 3 , kf or kf, 
and divide by twelve. One can write the dynamical zeta function as : 



CtW (1 - t)*i ' (1 - ^2 ) iV 2 • (1 _ t 3)N 3 (1 _ t r)W, - ■ ' (,{4 1 



The combination of the N r 's, inherited from the product (|34|), automatically takes into account the fact that the 
total number of fixed points of k r t can be obtained from fixed points of k\ , where p divides r, and from irreducible 
fixed points of k r e itself (see |36) for more details). A detailed analysis of this cycle decomposition (|34| ) for generic 
values of e will be detailed elsewhere [ p7[ . The previous exhaustive list of fix ed p oints (up to order twelve) can be 
revisited in this irreducible cycle decomposition point of view. The results of M yield : A x = 1 , N 2 = , N 3 = 
1 _ L A 4 = h_N 5 = 2 , N 6 = 2 , N 7 = 4 , iV 8 = 5 , A 9 = 8 , Aio = 11 , An = 18 . One actually verifies easily that 
( P9| ) and (p4| ) have the same expansion up to order twelve with these values of the A r 's. The next A r 's should be 
A 12 = 25 , A 13 = 40 , 7V 14 = 58 , A 15 = 90 , ■ • • 

It should be noticed that if one introduces some generating function for the real fixed points of k N , this generating 
function has the following expansion, up to order eleven, for e = .52 : 

Hl cal = t + t 2 + 4 1 3 + 5 1 4 + 1 1 1 5 + 4 1 6 + 15 1 7 + 13 t s + 40 1 9 + 31 t w + 67 1 11 + ■ ■ ■ (35) 

This series is a quite "checkered" one. Furthermore, its coefficients depend very much on parameter e . In contrast 
with generating function (pT|), the generating function iJ™ al has no universality property in e. This series does not 
take into account the topological invariance in complex projective space : it just tries to describe the dynamical 
system in the real space. This series i?J eal corresponds to the "complexity" as seen on the phase portrait of Fig. (||). 
One sees here the quite drastic opposition between the notions well-suited to describe transformations in complex 
projective spaces and the ones aiming at describing transformations in real variables. 



B. Dynamical zeta functions for a = 0, e non generic 

To further investigate the identification of these two notions (Arnold complexity-topological entropy), we now 
perform similar calculations (of fixed points and associated zeta dynamical functions) for e = 1 /m with m > 4 and 
e = (to — 1)/ (m + 3) with m > 7 odd. 

The calculations have been performed for e = 1/m for to = 4,5,7 and 9, giving the expansion of H e (t) up to 
order eleven. For to = 4 this gives : 

H 1/A {t) = t + t 2 + 4 1 3 + 5 t A + 11 1 5 + 10 1 6 + 22 1 7 + 29 1 8 + 49 1 9 + 71 i 10 + 111 t 11 + ■ ■ ■ (36) 

for m = 5 : 

Hi/ S (t) = t + t 2 + At 3 + 5t 4 + lit 5 + 16t 6 + 22t 7 + 37i 8 + 58 i 9 + 9H 10 + 144t n + ■ • • (37) 
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for m = 7 : 

H 1/7 (t) = £ + £ 2 +4£ 3 + 5£ 4 + ll£ 5 + 16£ 6 + 29£ 7 + 45£ 8 + 67£ 9 + lll£ 10 + 177£ n + --- (38) 
and for m = 9 : 

Hi/g(t) = £ + £ 2 +4£ 3 + 5£ 4 + ll£ 5 + 16£ 6 + 29£ 7 + 45£ 8 + 76£ 9 + 12H 10 + 188£ n + --- (39) 
All these expressions are compatible with this single expression of the C function : 

<l/m(t) = 1 _ t l S/ +t „ l+ 2 ( 4 °) 

We conjecture that this expression is exact at every order and for every value of m > 4. Again all the singularities of 
this expression coincide with those of generating function corresponding to the Arnold complexity (see Eq. (|l5|)). 

As far as functional relations relating £(£) and £(±l/£) are concerned, recalling (|30"|), one immediately verifies that 
£(f) corresponding to ([lO]) verifies the simple functional relation : 

- - t m+1 ■ Ci i (t) 

t^-C^)=C 1/m m, or: Cl/mm = tm+1 . (i/m ^Xl(t) + l (41) 

Actually (i/ m (t) has a very simple n-th root of unity form : 



Cl/mW — t ■ (1 - t m+1 ) 

Also note that when m is odd, and only in that case, Ci/m(t) a ^ so satisfies the functional relation 



(42) 



t rn+1 -(i /m (t) - -Ci /rn (-l/*) (43) 

No simple functional relation, similar to (^3|), can be deduced on H(t). 

Similar calculations can also be performed for the second set of non-generic values of e , namely e = (m— l)/(rn + 3) 
with m > 7, m odd. For m = 7, that is e = 3/5, one gets, up to order eleven, the same expansion as Eq. (pq) : 

ff 3/5 (t) = £ + t 2 + 4 £ 3 + 5 £ 4 + 11 1 5 + 16 i 6 + 29 1 7 + 45 1 8 + 67 £ 9 + 111 £ 10 + 177 1 11 H (44) 

suggesting, again, the dynamical zeta function : 

<8/.(i) - + (45) 

For m — 9, that is e = 2/3, one gets 

#2/3 (*) = £ + £ 2 + 4£ 3 + 5£ 4 + ll£ 5 + 16£ 6 + 29t 7 + 45£ 8 + 76£ 9 + 121£ 10 + 177£ 11 +••• 
A compatible zeta function could be[| : 

, , x 1 -£ 2 -£ n -£ 12 -£ 13 

- !_ t _ t2+t n (46) 

This form is not the same as Eq. (|4^), however it has the same poles. 

Comparing these rational expressions for the dynamical zeta function ((p9|), (p0|), ...), and the rational expressions 
for the generating functions of the Arnold complexity ((|l2|), jl 7 ^), ([Til), ...) for the generic, and non-generic, values 
of e , one sees that one actually has the same singularities in these two sets of generating functions (note that t 
has to be replaced by x 2 since k e is associate to transformation K 2 and not K). The identification between the 



7 The series is not large enough to confirm this form. A first simple analyzis seems to show that the next terms are ■ ■ ■ + 
296 1 12 +469t 13 +785£ 14 +•■■. 
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Arnold complexity and the (exponential of the) topological entropy is thus valid for generic values of e, and even for 
non-generic ones. 

It is worth noticing that, due to the topological character of the dynamical zeta function, these results are of course 
not specific of the y and z representation of the mapping (|3|) but are also valid for the (u , v) representation (Q) : 
in particular the exact expressions of the dynamical zeta functions (namely (29), (|40|)), remain unchanged and, of 
course, the denominators of the complexity generating functions are also the same for generic, or non-generic, values 
of e. 

The local area preserving property in the neighborhood of all the fixed points of fc™ previously noticed for a = 0, e 
generic, is also verified for these non generic values of e. 



C. Dynamical zeta functions for 

This (generic) identification is not restricted to a = . One can also consider mapping (Q) for arbitrary values of 
a and e and calculate the successive fixed points. Of course, as a consequence of the higher complexity of the a ^ 



situation (as shown in section II B, the complexity immediately jumps from 1.61803 • • • to 2.14789 • • • ) the number of 
successive fixed points is drastically increased and the calculations cannot be performed up to order eleven anymore. 
In the generic case, the expansion of the generating function H(t) of the number of fixed points can be obtained up 
to order seven : 

fff = 2i + 2t 2 + ll< 3 + 18i 4 + 47i 5 + 95< 6 + 212 i 7 + •■• (47) 

One has two fixed points for fc, no new fixed points for fc 2 , three sets of three new fixed points for fc 3 (giving 
3x3 + 2 = 11 fixed points), four sets of four new fixed points for fc 4 (giving 4x4 + 2 = 18 fixed points), nine 
sets of five new fixed points for fc 5 (giving 9 x 5 + 2 = 47 fixed points), fourteen sets of six new fixed points for fc 6 
(giving 14x5 + 3x3 + 2 = 95 fixed points). This expansion corresponds to the following order seven expansion 
for the dynamical zeta function : 

C(t) = 1 + 2t + 3t 2 + 7t 3 + 15i 4 + 32i 5 +69i 6 + 148i 7 + ••• (48) 

thus yielding to the following rational expression for the dynamical zeta function : 

C(t) P-^Hi±3 - g - ^ • (* + ^ with • t x 2 (49) 

This expression can also be written : 

C(i) . (' (50) 

l-t-(l+t) 2 

Let us recall the "alternative" zeta function (^). It verifies the simple functional relation : 

t 2 -c(t)-c(~t) = -C(-i/*)-Cf(V*) (5i) 

This new rational conjecture ( p9] ) corresponds to the following expression for Hit) : 

t-(2 + 3t 2 + t*) 

HAt) ~ (l-t 2 )-(l-t -2t 2 -t*) {b2) 

Comparing the denominators of Eq. ( [l9| ) and Eq. (|l6|), one sees that, like for the case a — 0, there is an identification 
between the Arnold complexity and the (exponential of the) topological entropy 

A = h (53) 

Hcuristically, this identification can be understood as follows. The components of k N , namely yjy and zjy, are of the 
form P N (y,z)/Q N (y,z) and R N (y, z)/S N (y, z) , where P N (y,z), Q N (y,z), R N (y,z) and S N (y,z) are polynomials 
of degree asymptotically growing like A . The Arnold complexity amounts to taking the intersection of the N-th 
iterate of a line (for instance a simple line like y = yo where yo is a constant) with another simple (fixed) line (for 
instance y = yo itself or any other simple line or any fixed algebraic curve). For instance, let us consider the N-th 
iterate of the y = yo line, which can be parameterized as : 
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Pn(vo,z) S N {y ,z) 

VN = 7T~7 \ ' ZN = 7M \ ' ( 54 ) 

Qn{Vo,z) T N (y Q ,z) 

with line y — ya itself. The number of intersections, which are the solutions of Pn(jjo, z)/Qjv(j/0; z ) — Vo> grows like 
the degree of P]y(yo,z) — QN{yo, z) • yo ■ asymptotically it grows like ~ X N . On the other hand the calculation of 
the topological entropy corresponds to the number of fixed points of k , that is to the number of intersection of the 
two curves: 

Pn{v,z) ~ Qn(v,z) -y = 0, Rn(v, z) - S N (y,z) ■ z = (55) 

which are two curves of degree growing asymptotically like ~ A^. The number of fixed points is obviously bounded 
by ~ A 2 N but one can figure out that it should (generically) grow like ~ A^ . This is fully confirmed by our exact 
calculations. 

The eulerian product Weyl-decomposition ( Jol] ) of the dynamical zeta function ( f49| ) corresponds to the following 
numbers of r-cycles : Ni = 2 , N 2 = , N 3 = 3 , N 4 = 4 , iV 5 = 9 , N e = 14 , iV 7 = 30 , N s = 54 , N 9 
107 , N w = 204 , Nn = 408 , N 12 = 25 , N 13 = 1593 , N 14 = 3162. 

D. Dynamical zeta functions for a/0 with e non-generic 

For a "non-generic" value of e when a^0, namely e = 1/2, the expansion of the generating function H(t) and of 
the dynamical zeta function read respectively : 

H? /2 (t) = 2£ + 2t 2 + lK 3 + 18t 4 + 47£ 5 + 95£ 6 + 198t 7 +••• 
Q /2 (t) = 1 + 2 1 + 3 1 2 + 7 1 3 + 15 + 32 1 5 + 69 1 6 + 146 i 7 +••■ 

A possible rational expression for the dynamical zeta function is for instance : 

r = i±*z£ = l+t-ci-t 6 ) 

l-t-t 2 -2t 3 -t 4 -2t 5 -t 6 -t 7 1 - t ■ (1 -t +t 2 ) ■ (1 +t +i 2 ) 2 V 1 

This last result has to be compared with (|l7|). 

For another "non-generic" value of e when a / 0, namely e = 1/3 the expansion of the generating function H(t) 
and of the dynamical zeta function read respectively : 

ff"/ 3 (t) = 2t + 2t 2 + lH 3 + 18f 4 + 42t 5 + 83t 6 + 177t 7 + --- 
C1/3W = 1 + 2 1 + 3 i 2 + 7 1 3 + 15 1 4 + 31 1 5 + 65 + 136 i 7 H 

A possible rational expression for the dynamical zeta function is for instance : 

c * (t) = l + t = 1+^ (57) 

^ 1/3W l-t-t 2 -2t 3 -t 4 -t 5 1 - t- (1 + t 2 ) ■ (1 +t +t 2 ) V 1 



This last result with has to be compared with (|18|). These results^ are again in agreement with an Arnold- 
complexity-topological-entropy identification. 

The local area preserving property in the neighborhood of all the fixed points of fc™ e previously noticed for a = 0, 
is also verified for a ^ for dTj) for generic values of e generic as well as these non generic values of e. 



8 However for the non-generic value of e , e = 3/5, we do not have enough coefficients in the expansion of the dynamical zeta 
function to actually compare it with Jig|). 
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IV. COMMENTS AND SPECULATIONS 



Based on analytical and semi-numerical calculations we have conjectured rational expressions with integer coeffi- 
cients for the generating functions of the complexity and for the dynamical zeta functions for various values of the 
parameters of a family of birational transformations. According to these conjectures, the Arnold complexity and the 
exponential of the topological entropy are algerbraic numbers. Moreover, these two numbers are equal for all the 
values of the parameters. 

From a general point of view, rational zeta dynamical functions (see for instance ]38|-[40[| ) are known in the literature 
through theorems where the dynamical systems are asked to be hyperbolic, or through combinatorial proofs using 
symbolic dynamics arising from Markov partition pJlj and even, far beyond these frameworks p2j, for the so-called 
"isolated expansive sets" (see |42| , [43|| for a definition of the isolated expansive sets). There also exists an explicit 
example of a rational zeta dynamical function but only in the case of an explicit linear dynamics on the torus R 2 /Z 2 , 
deduced from an SL(2, Z) matrix, namely the cat map |p|,|45f (diffeomorphisms of the torus) : 



A = 



2 1 
1 1 



B = 



1 
1 



det(l -z-B) 
det(l - z ■ A) 



{l-zf 



1 - 3 ■ z + z 2 



(58) 



Note that golden number singularities for complexity growth generating functions have already been encountered (see 
equation (7.28) in |Q or equation (5) in (4(|). In our examples, we are not in the context where the known general 
theorems can apply straightforwardly. The question of the demonstration of the rationality of zeta functions we have 
conjectured remains open. 

In the framework of a "diffeomorphisms of the torus" interpretation, the degree of the denominator of a rational 
dynamical zeta function gives a lower bound of the dimension g of this "hidden" torus C 9 jZ 9 . On expression ( |40| ) 
valid for a = and e = 1/m, one notes that dimension g grows linearly with m. The iteration of some birational 
transformations which density Abelian surfaces (resp. varieties) has been seen to correspond to polynomial growth 
of the calculations jy^]. Introducing well-suited variables 9i (i = 1, ■ ■ ■ g) to uniformize the Abelian varieties the 
iteration of these birational transformations just corresponds to a shift^ 9i — * 9i + n ■ r\%. For such polynomial 
growth situations, matrix A can be thought as the Jordan matrix associated with this translation, its characteristic 
polynomial yielding eigenvalues equal to 1. 

Many denominators of rational zeta functions encountered here are of the form : 1 — t-Y(t) where Y(t) is product 
of cyclotomic polynomials pO,pl| . We have encountered : 



Y(t) = (l + t) (resp. (l + t) 2 ) 



for 



l + t 



(i-ty 



i - 1 



(resp. a / 0) and 

for a ^ and e 



e generic , 
1 
2 



(resp.e = — ) 



More generally the rational dynamical zeta functions, or the rational functions G(q, x) encountered here, are of the 
form : (1 + X{z))/(1 - Y(z)) (for G) or (1 - X(z))/(1 - Y(z)) (for Q where X{z) and Y{z) have some kind of 
decomposition on cyclotomic polynomials : 



X(z) 



z r ■ H m P^\z) with P£'(z) cyclotomic polynomials 



(59) 



This is particularly obvious on expressions ( J14| ) but also on expressions (p. 3D, or (|56|), or even (|40|). We do not have 
yet any Z-adic cohomology interpretation (see for instance (3(| page 453) of this cyclotomic polynomials "encoding" of 
the zeta functions or the complexity functions G(q, x) . Most of these rational expressions for zeta functions satisfy 
very simple functional relations but one also expects, for ( |56|) or (|5^ ) for instance, more involved but, still simple, 
functional relations similar to the ones obtained by Voros in |52]. Many of the generating functions G(q, x) can also 
be seen to satisfy simple functional relations relating G(q, x) and G(q, 1/x). This will be detailed elsewhere^]. 

In practice it is numerically easier to get the generating functions of Arnold complexity than getting the dynamical 
zeta functions. If one assumes the rationality of the dynamical zeta function and the identification between Arnold 
complexity and (exponential of the) topological entropy, getting the generating functions of Arnold complexity is a 
simpler way to "guess" the denominator of the dynamical zeta functions. 



This "diffeomorphisms of the torus" interpretation is quite obvious on figure 2 of 
10 For instance the generating function of the degrees g(x) given by equation (5) in 



46] 



verifies g(x) + g(l/x) 
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The analysis developed here can be applied to a very large set of birational transformations of an arbitrary number 
of variables, always leading rational generating functions f44|j4{| . Moreover, these generating functions are always 
simple rational expressions with integer coefficients (thus yielding algebraic numbers for the Arnold complexity). 
They even have the previously mentioned "cyclotomic encoding". At this point, the question can be raised^ to see 
if the iteration of any birational transformation of an arbitrary number of variables always yields rational generating 
functions for the Arnold complexity. We have even found rational generating functions of Arnold complexity for 
rational transformations which are not birational (see (7.7) and (7.28) in pi|] ) : any proof of these rationalities should 
not depend too heavily on a naive reversibility of the mapping |p3| . 

We have also calculated Lyapunov exponents pm in order to study the metric entropy. These numerical calculations 
will be detailed elsewhere [[34| for transformation^ for e = .52. These results share quite small values of the Lyapunov 
exponents, the largest of which being much smaller than the topological entropy. We thus infer that, in this very 
example, the metric entropy is much smaller than the topological entropy. We have here an opposition between 
topological concepts originating from complex projective spaces and the metric concepts of real analysis. The "non- 
topological" complexity measures do not seem to be able to identify with the previous topological and algebraic 
quantities. On the birational examples studied here, the metric entropy does not seem to share the same algebraic 
values as the topological complexity measures. 
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APPENDIX A: CORRESPONDENCE BETWEEN TRANSFORMATION K ACTING ON Q x Q 
MATRICES AND THE (U , V) TRANSFORMATION K a , e 



In previous papers p0|j2l| ] , it has been noticed that the successive even iterates of K , acting on an initial matrix 
M, actually belong to a plane which contains matrices M , K 2 (M) and K A (M), or equivalently a fixed matrix P 
(see (§) : 

K 2n {M) = cq ■ M + ci • K 2 {M) + c 2 • P 

for any integer n (even relative integer). In fact one even has the following property : any point of the plane containing 
M, K 2 (M) and P is transformed, by the even iterates K 2n , into another point of this plane. This can be used to 
define the two dimensional mapping k e>a compatible with mapping K . Let 's M be an arbitrary q x q matrix, and 
let us define^ : 

u n (M) = x 2 ■ x± ■ x 6 ■ ■ -x 2 n-2 ■ Ui(M) , v n (M) = x 2 ■ x 3 ■ x± ■ x 5 ■ x 6 ■ ■ ■x 2n -i ■ vi(M) (Al) 
where: Ul (M) = ^, « 1 (M)=-^, and: a(M) = , e(M) = (A2) 

P2 -^2 \ 2 \ 2 

where x n = det(K m (M)) ■ det{K n+1 (M)) with K = t ■ I, and T(M) = M^ 1 and[] : 

x x 2 ■ (xix 3 - 1) xi ■ (x x 2 - 1) 

^i — --, j * 2 — , 

x 2 — 1 x% — 1 

X\X 2 Xz + X\X 2 — X\ — 1 XqX\X 2 + X Xi ~ Xo ~ 1 

Pi = ; j P2 = : 

x 2 — 1 Xl — 1 
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in equation (6.35) u„/v n should be replaced by v n /u n yielding (see (A2) and (A2)) : v n /u n 



1 After 0. 
2 Note a miss-print in 

X- A ■ -X;, ■ X 7 ■ ■ ■ X 2n -1 ■ V\/u± . 

13 Note that the Ai's and piS are not exactly the same as the ones given in in equations (6.13) and (6.14) : the A^'s and 
Pi's in j2(| are homogeneous expressions and the A^'s and p^'s we introduce here are inhomogeneneous true invariants which 
can be deduced from the ones in dividing them be qo or qi . 
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Then one has the K 2 invariance of a and e : 

a{M) = a{K 2 {M)) , e(M) = e(K 2 (M)) (A3) 
and fc Q . £ can be seen as a representation of K 2 : 

(u(M),v(M)) = k a>e (u(K 2 (M)),v(K 2 (M)) (A4) 
where a and e are precisely the values given by Eq. ( |A2| ). Transformation /co,, e reads : 

k ae : (u,v) — > (U , V) = ( , U ■ a + 1 + e ) (A5) 

11 u 

In the a = case, this transformation simplifies and one can introduce new variables y and z given by : 

y = v - 1, z = — — —, - ^ (A6) 

u ■ (v — 1) 

With these new variables, k a . e reads : 

k e : (y,z) — » (z + l-e, 7) (A7) 

2+1 

For a = e = 0, transformation fc Q £ is integrable |2(J the invariant being (see (6.38) in p0[| ) : 

I = (A8) 
u 

This algebraic expression is of course only well-suited for e = 0. The variable z amounts to considering X/ (v — 1) for 
arbitrary e's. 

APPENDIX B: FACTORIZATION SCHEME FOR (§) 

For q x q matrices (q > 3) the factorizations corresponding to the iterations of K read: 

A--**), h-^jgL, M 2 = S*>, /3 = Mp>. M S = ^>, 

7i / 1 71 ' J 2 J 2 

= rfet(M 3 ) = ff(M 3 ) = rfef(M 4 ) _ tf(M 4 ) 



f9-l f f9-2 ' f g-2 f q-3 ' " f2 f ?-l r f q-2 ' ^ f f q-2 f q-3 ' 

J 1 " J2 ■ J 3 7l ■ 7 3 h J2 " 73 J4 71 ' J2 ' 74 



drf(M 5 ) _ 7^(M 5 ) . _ det(M t 



Mi) 



76 — „ g _2 f 2 f f9-2 ' 6 _ f g-3 /• rq-2 rq-3 ' f fl^ 2 $2 fl^ 1 f fl^ 2 

7i • h ' 7 3 ' 74 • 7 5 7i ' 72 • 7 3 • 7 5 71 ' 72 ' 7 3 ' h ' 75 • J 6 



M 7 = - ? ... (Bl) 



#9— J f f9-2 f8 
J2 ' 73 ' J4 ' Jq 

and for arbitrary n : 

det(M n ) = f n+1 ■ (/r 2 • fn-l ■ ft-l ■ fl-3) ■ tft-l ■ fn-B ' ft-\ ' /^- 7 ) ' ' ' /x" , (B2) 

7Y(A7„) = M„ +1 ■ (/r 3 • /r' • /n-s) • (/^4 • /re • fn-i) ■ ■ ■ /f" (B3) 

where fi n = q — 3 for n = 1 (mod 4), /i„ = for n = 2 (mod 4), /i„ = q — 2 for n = 3 (mod 4) and \x n = 1 for n = 
(mod 4) and 8 n also depends on the truncation. Factorization relations independent of q, occur : 

K(M n ) = A^+i 

det(M n ) fn+l • fn ' /n-1 ' fn-2 ' fn-3 ' fn-4 ' ' ' 

Let us introduce p0| , ^l[ the generating functions a(x) and (3(x) of the degree of the det(M n )'s and /Vs. Their 
exact expressions read: 

()- 9 + 9 2 ■ a; • (1 + s 2 ) g-x.(l + a; 2 ) 

1+x (1 — x)(l + x)(l — x — x 3 ) ' 1 — x — X 3 

It is clear that one has an exponential growth of exponents a„'s, /3„'s, /i n 's and f„'s: these coefficients grow like A™ 
where A ~ 1.465 • ■ • 

This displays the "generic" factorization scheme. However, on various subvarieties like codimension one subvariety 
a = , the factorization scheme can be modified as a consequence of additional factorizations occurring at each 
iteration step, thus yielding a smaller value for the complexity A . 
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1. Factorization scheme for a = 0, e generic 



For a — the previous factorization scheme becomes for 3x3 matrices^ : 

f j t(M \ m v t at \ f detjMt) det{M 2 ) K(M 2 ) 
ji = det(M ) , Mi = K(M ), f 2 = -, , M 2 = K(M 1 ), j 3 = — — — , M 3 = — 



deUMs) M tc(M \ f det{Ah) M K{M4) 



f f f i * \ *5 / ' JO f2 (2 f2 f ' f f f ' 

Jl ' J2 ■ J 3 Jl ■ J2 ' 73 ' 74 Jl ' /2 ' /3 

_ def(M 5 ) _ K(M 5 ) _ det(M 6 ) K(M 6 ) 

J 6 f „2 f f f ' -™6 /• ? 77 f2 f f2 f2 f2 f ' 7 f f f f 

Jl ■ H 73 " 74 ■ J5 72 Jl ■ J2 ■ J 3 ■ Ji ■ J 5 ■ J6 Jl • /3 ' J4 • /5 

rfef(M 7 ) Jf(M 7 ) dei(M 8 ) 

78 — t j j 72 — 1 1 T > Ms = 7 ' J 9 



f f f f2 f f f 1 ° f ' J- J f2 f2 f2 f f2 f2 f2 f ' 

Jl ■ 72 ■ 73 • J 4 ' 75 ' 76 ' 77 74 7i • 72 ' J3 ' 74 • 75 " 76 " 77 ' 78 

g(Mg) det(M 9 ) 

9 f f t f t f ' JW ff2ffff2fff' \ DU ) 

7l ' 72 ■ 73 • 75 ■ 76 • 77 71 ' 72 • 73 • 74 • 75 • 76 • 77 • 7s • 79 

and for arbitrary n : 

det(M„) = / n+1 • (/„ • /2_ x • / 2 _ 2 • fl_ 3 ) ■ (/„_ 4 ■ / 2 _ 5 • / 2 _ 6 • /»_ 7 ) • • • (B7) 

K(M n ) = M n+1 ■ (/„_! • / n _ 2 • /„_ 3 ) • (/„_b • f n - 6 ■ fn-7) • • • (B8) 

for n even and : 

det(M„) = f n + 1 • (/„ • /„_! • / n _ 2 • /*_ 3 ) ' (fn-4 ■ fn-6 • /»-6 ' /'- 7 ) ' • ' (B9) 

tf(M n ) = M n+1 • / n _ 3 • /n-7 ' /n-U ' /n-15 ■ /n-19 ■ ' ' (BIO) 

for n odd. 

The exact expressions of the generating functions a(x) and (3(x) rcadp'l : 

a(x) = + , where: 0(x) = 3 . * ' C 1 + * + * 3 ) = _ 3 + 3 ■ (1 + x)/(l - .t 2 ~ x 4 ) (Bll) 

1 + a; 1 — x z I — x z — x^ 



It is importa nt t o note that factorization scheme ( JBq ) is actually stable, but of a slightly more general form, as 
compared to (Bl), or the ones described in Q] : recalling the generating functions n(x) and (f>(x) of the exponents 
that occur in the factorization scheme (see equation (8.6) and (8.10) in |t4|), one must now introduce two sets of such 
exponents generating functions, rji, <f>\, 772, <p2, in order to keep track of the parity of n, and even split these four 
functions into their odd and even parts : 

T712 = {vi(x)+Vi{-x))/2, mi = fai(aO-T?i(-a:))/2, (B12) 
V22 = (V2(x) +t] 2 {-x))/2, mi = -rj 2 (-x))/2, 12 = ... 

We must also decompose a(x) and f3(x) in odd and even parts: 

a(x) - a(-x) a(x) + a(-x) /3(x) - f3(-x) (3{x) + /3(-x) 

"lW = o > Ot2[X) = , Pi(x) = , P 2 {X) = 



namely 



3 ■ x* ■ x 2 + 1 3-x 

3- (l + 2x 2 + 2a; 4 ) , , 3-a;- (2 + x 2 + a; 4 ) 

a 2( x ) = 7i tVTI 2^ IT' a i( x ) 



(1 - x 2 ) (1 - x 2 - x 4 ) ' 1 w (1 - x 2 ) (1 - x 2 - x 4 ) ' 



14 These r esult s can straightforwardly be generalized to q x q matrices, they are just a bit more involved. 
15 Result (Bll) corresponds to a very simple expression for another generating function introduced in fl4| ], namely the function 



p(x) (see for instance equation (8.12) in Q) 
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Instead of functional relations (8.6) and (8.10) in Q], one now has the following relations 



ai(x) - 2 • x ■ a 2 (x) + 3 • x ■ (t/i2(x) ■ (3 2 (x) + 7/ii(x) ■ /?i(x)) = , 
a 2 (x) — 2 • x ■ eti(x) - 3 + 3 ■ x ■ (7722(2;) ■ Pi{x) + 7/21(1) • (fa(x)) = , 
x ■ ai(x) - (3 2 {x) - (021 (x) ■ Pi(x) + <p 22 (x) ■ (3 2 {x)) = 0, 

x ■ a 2 {x) - (3x{x) - {4> u (x) ■ fo{x) + <pi 2 {x) ■ /?i(x)) = (B13) 
where the odd and even part of the exponents generating functions T)i(x), (/>i(x), ^2(^)7 4>2 (#), read : 

7/12(2:) = -, T, 7/n (x) = J' ^22(2;) = 0, 7/21 (x) 



1 -x 4 ' ' v ' 1 -x 2 ' ' v ' ' ' v ' 1 -X 4 ' 

1 / x x • (2 x 2 + 1) . . n x 2 , . , x , , , x 2 • (2 x 2 + 1) 

M*) = — — . M*) = 2 7^2' ^iW = ^3^' M*) = — f3^4 — . 

Period four in the factorization scheme ( p36| ) corresponds to the occurrence of a 1 — x 4 = singularity for these 
exponents generating functions. 

The "stability" of factorization scheme ( pl| ) corresponds to the following (n — > n + l)-property : the exponents 
of the /„'s occurring at the m-th step of iteration are also the one's at (m + l)-th step of iteration, the /„'s being 
changed into f n +i '■ at each new iteration step one only needs to find the exponent of fx (if any). The "stability" 
of factorization scheme (B6) is a straight generalization mod. 2. of the previous property : the exponents of the / n 's 
occurring at the m-th step of iteration are also the one's at (m + 2)-th step of iteration the / n 's being changed into 

fn+2- 



2. Factorization scheme for a ^ 0, e non generic 

Let us come back to a / with the non-generic value e = 1/2. We consider here a = 396/6095 ~ .06497128. 
The factorization scheme reads : 

h = det(Mo) , Mi = K(M ) , f 2 = de * (Ml) , M 2 = if (Mi) , / 3 = <^MA , m 3 = K(M 2 ) , 

71 Jl ■ 72 

_ det(M 3 ) _ K(M 3 ) _ det(M 4 ) _ K{M A ) 

74 — ~F> 7 7~ ' 1V1a ~ 7 I 75 — To 72 7 7~ ; M 5 — ~? 7 , 



£% £ f ' 4 -f ' '° f2 J-2 J- J- ' J- J- 

Jl • J2 ■ 73 Jl Jl ■ J2 " 73 ■ 74 Jl ■ J: 



2 



_ rfe*(M 5 ) „ _ JT(M 5 ) _ rfet(M 6 ) _ jf(M 6 ) 

,/6 f f 2 (2 f f i M 6 f f i J 7 f f f2 f2 f f ' M7 f f 

J 1 ' J 2 ' J 3 ' 74 ' J5 /2 ' /3 Jl ' J2 ■ J 3 ■ J4 ■ J5 ■ J 6 J 3 ' J 4 

_ det(M 7 ) K(M 7 ) det(M s ) 

J8 — -p? — 7 7 7? — 79 — 7 7" ! M 8 = ~, 7 7" , J9 



/l • /*2 • /3 • /f • /| • /e • h ' /l " /l ' /5 ' /i • /| • /3 • /4 • /| • /I ' h ■ /s 

_ g(M 8 ) _ rfet(M 9 ) *T(M 9 ) 

M 9 - 1 7 7 7" > -ho - 72 — 72—7 7 fl — 72 — 7 T ' M W - 



fx ■ /2 • h ■ /6 ' /i • /I • /I • fa ' /s • /I • /I ' /s ■ h ' h ■ h ■ /e • /7 

. det(Mio) ^(M i8 ) 

/ll - ~, 7 72 — 72 — 7 7 72 — 72 — 7 7" > ' ' ' M 19 



/20 — 
M 20 - 



/l ' /2 - /| - /I ' /*5 ' 76 - /| " /| ' /9 - /lO ' $3 ' f§ ' fl ■ fg, ■ 111 ' fl2 ' fl5 ' /l6 

de7j(Mi 9 ) 

f\ ' h' h ' h ' fe ' f$ ' fao ■ fll ■ fl2 ■ fl3 • fu ■ fl5 ■ fie ■ fl7 • hs ■ fl9 ' 

K(M 19 ) 

fl ■ f5 ■ fg ■ fa ■ fl2 ■ fl3 ■ fa 6 ■ fl7 ' 



f _ det(M 20 ) 

J 21 f2 f f2 f j-2 f2 f2 f f f2 f2 f f f2 f2 f f 1 V ^) 

Jl ■ J3 ■ J 5 ■ J7 ■ Is ' fa - JlO ' hi ■ J 12 ■ J 13 ' 7l4 ' 715 • 716 ' J 17 " 7l8 ' 719 ' 720 

Up to the thirteenth iteration one has the previously described ( n — > n + l)-property, but this property is broken 
with /i 5 in favor of the (n — > n + 2)-property encountered with (B6). The previously introduced odd-even-parity 
dependent exponents generating functions rjij (x) and (x) now read : 

™15 

10 1 „12 „ /„\ „3 , J7 , „11 



7/12 (x) = X 2 + X b + X iU + X 12 , 7/n(x) = X- 5 + X 7 + X 11 + 



1-x 4 ' 
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7722(2; 

011 (x 

012 (a; 

021 (H 
022(52 



X- 

(1 



x 6 + x 10 



2x 6 + 2x' 



,.14 



2a; 2 



,.14 



1-X 4 

x + 2x 3 + 2x 7 

„4 1 _6 



h 55 

2 x s 



2x xi 

2 a; 4 



2x 



11 



7721(3;) = x 3 + x 7 + X 11 
x 5 + 2x 13 , 
x 6 + 2x 8 +x 10 - 
js , (1 



2a; 



12 



-2i 2 ).x 13 
1 - a; 4 



■2a; 1 



a; + 2 x + x" 

from which one deduces, from relations flB13| ), the rational expressions of the ctj's and /3j's : 

3 -a; 2 ■ (1 + a; 2 ) 



/3 2 (x) 
a 2 (x) 

cti(x) = 3 • x ■ 



(1 -x 2 ) • (1 



2 a; 6 



2x x 



x L 



3- (1 + a; 2 ) • (1 + a; 4 ) ■ (1 + x 8 ) ■ x 



= 3- 



1 - a; 2 - a; 4 - 2 x 6 
1 + 2 a; 2 + 5 a; 4 4 



- x° 
4 a; 6 



2 a; 10 
■5 a; 8 - 



£.12 

4x 



X" — x 
10 



14 ' 

5 a; 12 



5 a; 1 



3 x 
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.8 



(1 - a; 2 ) • (1 - a; 2 - a; 4 - 2 a; 6 - x 
(2 + 4 a; 2 + 4 a; 4 + 5 x 6 + 4 a; 8 + 5 x 



2x 
,,10 



in 



4 a; 1 



-a; 14 ) 
4x 14 ) 



2 a; 6 



(1 - x 2 ) • (1 -x 2 -3 
yielding the rational expression for P(x) : 

3 • x ■ (1 + x + x 3 - x 16 ) 



2a; 



10 



,•12 



a; 14 ) 



p(x) 



1 - 2 a; 2 



,.!(> 



(B15) 



3 ■ 



X ■ (1 + x 2 ) ■ (1 + x - x 2 + . 



,.10 



,.12 



„14\ 



(1 - x 2 ) ■ (1 - x 2 - x 4 - 2 x 6 - x 8 
The complexity growth corresponds to the (smallest) root of : 

1 - a; 2 -x 4 -2x 6 - x 8 -2x 10 - x 12 - x 14 



2x 10 -x 12 ~x 14 ) 



= 



(B16) 



These results have also been checked using the previously depicted semi-numerical complexity growth evaluation 
method for e = 1/3 and a = 396/6095 ~ .06497 ■■• . The following value for the complexity has been obtained : 
A ~ 1.46199 in good agreement with the exact algebraic value deduced from (B16), namely : A ~ 1.46188 • • • (to 
be compared with the gene ric algebraic value of A , A ~ 1.4655 • • • associated with 1 — x — x 3 =0). 

The singularities of (|B15J) are in agreement with the dynamical zeta function calculated for these values of a and e: 

c(t) = 1 + _ l+f-(l-t 6 ) 



1 - t - t 2 - 2t 3 - t 4 - 2t 5 - t e ~ t 7 



1 



t-(l -t +t 2 )-{l +t +t 2 f 



(B17) 



These calculations can also be performed, for a 7^ 0, for the other non-generic value of e : e — 1/3. As far as 
the factorization scheme is concerned one gets exactly the same scenario as the one described in (B14), the breaking 
of the (n — > n + l)-property and the occurrence of a (n n + 2)-property taking place with /n instead of /15 
previously. For e = 1/3 and, for instance, for a = 237/6095 ~ .038884 • ■ 

/3(x): 



p(x) 



3-x-(l + x + x 3 — x 



12-1 



l-2x 2 



X 6 + x 8 



,.12 



one gets the following expression for 

(B18) 



3-x-(l +x 2 )-(l +x 



x 10 ) 



(1 



(1 



2x 6 



x 10 ) 



Again these results have been compared with the complexity growth deduced from the semi-numerical method, for 
e = 1/3 and a = 237/6095 ~ .038884 • • •. We have obtained the following value for the complexity : A ~ 1.44865 
in good agreement wit h the exact algebraic value deduced from ( B18| ), namely : A ~ 1.44717 • • •. 

The singularities of ( B1S| ) are in agreement with the dynamical zeta function calculated for these values of a and e: 



C(*) 



1 + t 



1 + t 



1 - t - t 2 - 2 1 3 - i 4 - 1 5 



1 -t-(l + < 2 )-(l +t +t 2 ) 



(B19) 
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APPENDIX C: CHOICE OF A INITIAL MATRIX CORRESPONDING TO GIVEN VALUES OF e AND a 



We present, in this section, a possible choice of an initial 3x3 matrix corresponding to a prescripted value of a 
and e. From the results of Appendix |A|, one has : 



a 



and : 



(x 3 X\ X 2 + X 2 Xl — X\ — 1) ■ {x 2 X Xi + Xi x - x - 1) 
(x 2 - 1) • (x x 2 - 1) -xi 

_ (xi x 3 - 1) ■ (1 - Xl) ■ Xq ■ x 2 

(1 - X 2 ) ■ (X X 2 - 1) • Xl 



1 



(CI) 



(C2) 



In order to perform our complexity growth calculations to get the factorization scheme of the transformation, one 
needs to iterate a non-trivial, initial matrix as simple as possible, in the a = case and for specific values of e 
(e = .52 , e = 1/m , ...). Actually, let us consider a matrix of the form : 



1 3 x 
5 2 y 
-4 8 z 



(C3) 



The a = condition factorizes as follows : 

_ (y + 5) (x + 3 + (2x + 61 - 11 y + 2 (a; - y - z) 
(y-5) 2 (z + 4:x)(2x- lly-61 + 2z) 

On the other hand, expression of 1 + e is also very simple since it also factorizes : 



1 



(x — z — 5) (5 x + 5 z + 3 y) 
(y - 5) 2 (z + 4 2) (2 x - 11 y - 61 + 2 z) 



(C4) 



(C5) 



APPENDIX D: THE POLYNOMIAL TO FIND THE FIXED POINTS OF Kl 

The fixed points of k^ 1 can be found as suitable pair of roots of two polynomials P(z) and Q(y). The number of 
pairs of roots being relavitely small (degre(P) x degre(Q)), it is straightforward to check which are the admissible 
pairs. For e = 13/25 and N = 9, the two polynomials happen to verify P(x) = Q(—x), We give below the expression 
ofP(z): 

P{z) = 314414322376251220703125 z 18 + 1358269872665405273437500 z 17 

+75268905252456665039062500 z 16 + 281939167586425781250000000 z 15 
+4712354272080487976074218750 z 14 + 14702451771291308349609375000 z 13 
+115459295503780457067138671875 z 12 + 289162068299094274224609375000 z 11 
+235039074495145372852311328125 z w - 28423190864054603531819812500 z 9 
-129391896463704494904550698750 z s - 47468841855664870004702580000 z 7 
+9768520701929861757756144700 z 6 + 8841684508557014424153308400 z 5 
+1497468490621088327339020023 z 4 + 77417791834794939443209320 a; 3 
+14196266775922682562956676 z 2 - 525991376147246600507280 z 

+4602174329226460987728 (Dl) 

The actual value of z for the fixed point ony + z = 0is: z~ -.4956845 + .003449852 ■ J. Polynomials P(z) 
and their partners Q(y) corresponding to the fixed points of fc 10 , fc 11 , and fc 12 , are available in |]34]] as well as their 
respective pairings of roots. 
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